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Abstract 

Iterative solution of QED evolution equations for non-singlet electron structure func- 
tions is considered. Analytical expressions in the fourth and fifth orders are presented in 
terms of splitting functions. Relation to the existing exponentiated solution is discussed. 
PACS: 12.20.-m Quantum electrodynamics, 12.20.Ds Specific calculations 

§ 1 Introduction 
a^ 

0> ; In this paper we are going to discuss properties of the QED non-singlet splitting functions. 
\ Some details of derivation of the fourth and fifth order approximations are given. The related 
subjects concerning the QED structure functions (SF) themselves are touched very briefly, 
mainly to show, how do the higher order splitting functions enter into the SF. An extended 
^ , discussion about the SF can be found in papers [|l], ||, H, ^ and references therein. In Ref. 

a fifth order perturbative solution for the non-singlet SF was derived within the ad hoc expo- 
I nentiation procedure 0. 

■ The Dokshitzer-Gribov-Lipatov-Altarelli-Parisi evolution equation for the non-singlet 

electron structure function reads 

V^'iz, Q^) = 5(1 - + / ^ ^ / - P«(x)P^^(^, q% (1) 

J ITX q"^ J X X 

where m is the electron mass; P*^^^ is the first order non-singlet splitting function; a{q^) is the 
QED running coupling constant. Here we are going to consider only the electron contribution 
to vacuum polarization: 

a(g2) = ^ (2) 

1 - ^ In 4 

It is worth noting that only the one-loop approximation (re-summed) gives the leading log 
contribution, while higher orders provide only next-to-leading corrections. The account of the 
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running coupling constant in Eq. (|1|) is associated with the leading log radiative corrections due 
to pair production in the so-called non-singlet channel. Usually pair production is considered 
separately from the pure photonic correction, because of different conditions of registration. 
Nevertheless, we are going to evaluate the general equation, and then point out the part, 
corresponding to pair production. An account of other non-singlet pair contributions due 
to muons, r-leptons, and hadrons can be done within a certain approximation in the final 
formulae. 

The singlet pair production mechanism in the leading logarithmic approximation is de- 
scribed by the so-called singlet electron structure function (see [||, ||, §]); we are not going 
to discuss it here. 



1.1 Splitting function 

Function 

p(i) 

(z) in Eq. ([1|) is the so-called splitting function. It serves as a kernel function 
in the evolution equation. The function (multiplied by the proper coefficient) gives the first 
order leading log correction to the probability, which is just P^'^(z, Q^), to find an electron 
(quark) with energy fraction z in the initial electron (quark): 

P<..(.) . P^W ^ (i±£)^ - ^(1 - .) /d. i±^. (3) 



A simple derivation of the function can be found in Ref. [0. 

The splitting function is a generalised mathematical function. For practical applications we 
prefer to use the following definition, which allows to avoid explicitly the mutual cancellation 
of infinite quantities during numerical computations: 



P«(z) = lim{2(lnA + ^)5(l-^) + i±^e(l-z-A)} 

= P^^^5il-z) + P^'^eil-z-A). (4) 

As could be seen from Eqs. (^,p, the function satisfies the following normalisation condition: 

1 

JdzP^^\z) = 0. (5) 



This condition is a manifestation of the Kinoshita-Lee-Nauenberg theorem P]: it provides 
the cancellation of mass singularities. 



2 Iterative solution 

In QED the evolution equation can be solved to any desired order of perturbation theory by 
means of iteration. The initial approximation for the structure function is just 6{1 — z). The 
first iteration gives 

P^^(^, Q') = 6(1 -z) + ^P«(^) + Oia'), P=-iL-l), L = In ^ , (6) 
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where L is the so-called large logarithm. 

On the next step of the procedure we need to calculate the integral 



dx 



X 



(7) 



This is a typical Mellin convolution 



P("+l)(^) = / ^P^^\x)P^'"\^) = j^^^j P^^\XI)P'^''\X2)S{Z - X1X2). (8) 



In this way step by step we get the solution for the evolution equation to the fifth order: 

V^'{z,Q') = V^'iz,Q')+V^,Uz,Q% (9) 

(10) 



n=l "'■ 



V. 



NS 



+ -P(i) 
27 



4 



3 27 



+ 



6 54 



4 



+ 



16 



^p(4) I ]_ p(3) , 1^^ p(2) , 

18 54 81 405 



p(i) 



0{a'). (11) 



We denoted by index 7 the pure photonic part of the SF. The other part describes pair 
corrections, and, starting from the third order, with possible simultaneous photon radiation. 
Recently we considered the numerical impact of the higher order pair corrections to electron- 
positron annihilation in paper 

At the level of the non-singlet structure function the condition (^) reads 



_L 1 

j dz p^^(z, = 1, y dz pW(2) = 0, 



n= 1,2, 



fl2) 



This condition has also a trivial probabilistic meaning: the sum of probabilities of all allowed 
emission processes is unit. In other words, one is always able to find an electron in the initial 
electron. 

Now we return to the splitting function properties. Using prescription (^) we can represent 
the integral of the product of our two generalised functions in the form 



p(»+i)( 



Pi"+')5(1 -z)+ p^"+^^e(i - ^ - A), A ^ 0, 

:,(n)jZs 



e 



X 



PS\z)Pi^^ + P^^P^-\z)+ j ^^P^\x)P, 

J X 

z/{l-A) 

The A-part of the splitting function can be obtained from the condition (|5|JI^): 

l-A 



P 



(n+l) 



d^P^ 



Z). 



(13) 



(14) 
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Instead of the above trick we can use, as was discussed in |T[], the known solution of the 
evolution equation in the soft limit [10(1: 



p (1 - z)"/^-^ fpf3 

2 r(i + /5/2) 



(15) 



where C is the Euler constant (C ~ 0.57721566). In order to obtain the A-part of a splitting 
function of the desired order we have to integrate Eq. (^) over the interval 1 — A < ^ < 1, 
and then expand into a series in a: 



1-A ^ ' ^ 



(16) 



Now we have to expand the exponent and use also the following formula (see Appendix): 

'C(4) 



T{l + /3/2) " 2V2 



+ 



5 V2 



C(5) 



3 V2 
6 V2 



16 V 2 

C(2)C(3) + 0(/3^). 



(17) 



The second and third order splitting functions are well known (see paper |l| and references 
therein) and used in many various applications. For the sake of completeness we put here the 
expressions: 



,(2) 



^^^^ ^21n(l-^) -ln^ + ^ 



1-z 



4 InA 



1 + z 



\nz — 1 + z 



24 



1+z^/l 



4C(2), 



In'il-z) 



;i8) 



3 1 1 
- ln(l — z) In 2; ln(l — z) -\ In^ z 

4 ^ '2 ^ '12 



I- z\2 

^ Inz + - -C(2) ) + 6(1 + z) ln2ln(l - z) - 12(1 - z) ln(l - z) 



8 
3 



32 



+ -(5-3z)ln2-3(l-^) 



-(1 + z) z + Q{1 + z)U2{l - z) 



P. 



(3) 



In A + - 24C(2)fln A + ^ U 16C(3). 



4 



4 



(19) 



The definition of the Riemann (^-functions, dilogarithm, and other special functions are given 
in Appendix. 

By means of the convolution procedure (0,0) we found 



144 
1 



rl+^ 

\l-z 



- \n\l -z) + - -z) + {-- -C(2) ln(l - z 



9 



- ln^(l — z)hiz -\ — ln(l — z) Ir? z ln(l — 2;) In 2 

3 9 2 



In^ z In^ z 



108 



+ 



12 



+ 



V 3 16 



1 2 3 C(2'\ 4 

\nz-- ln.Li,(l -z)- -81,2(1 " ^) + ^ " ^ + gClS) 



4 



__ln2(l_^)__ 
l + z 



6 



ln(l -^) + 



1 + ^,2/-, \ 1 

— ^ — 111 (1 — z)mz 



1 /-. \i 2 5 — 32; , 7(1 + 2;) , o 52; — 11 , 2 
ln(l - 2;) In^ 2; H — — ln(l - 2;) In 2; + ^ In'^ z H r^— In^ z 



12 



^64 



144 



+ 
+ 



12 

43 — 52; 1 + 1 + -2 , x-r . X 1 + 2;^ . ,^ , 
C(2) In^ + -—- ln(l - ;2)Li2(l - z) ^^3(1 - z) 



288 

l + z 



6 



6 



Si,2(l -z) + -^Li2(l -z) + ^^C(2) ^ 



12 



144 



(20) 



p(4) 



16(^ln A + - 96C(2) (^In A + + 128C(3) (^In A + ^) + 24C(4). (21) 



At the next step we got 



P§\z) 



+ 
+ 



+ 
+ 



+ 



I - z 



UAi-z) + \\n\i-z) + ^^ ' 



16 



C(2)) ln2(l - z) + ( 
- C(2) + ^C(3)) ln(l -z)-hnz \n'{l - z) + z -\ In z^ \n\l - z) 



-1 \^'z + \ In^ z - (^Li2(l -z) + \- In z - ^S,,2(l - z) ) ln(l - z) 



— In^ z - — In^ ^ + f — - ic(2) + — Li2(l - z)\ In^ ^ + (-— + ^((2) 
)80 72 V16 9^^ ^ 18 ^ 7 V 32 2^^ ' 

C(3) + lSi,2(l - 2;) - lLi2(l -z) + ^Li3(l - 2;)) In z + h.il{l - z) 

(1 + ^) 



10 
4 

9 

^^^'^(^-"^ + 256 -8 



^C(2) + ^C(3) + ic(4) 



9 



In 2; 



2(1-2;) 



9 

l-z' 



ln3(l - z) + 



12 



5-32; (1+^ . 

In2; + ^ — - — ^Li2(l — z) 



12 



\ 1 2/-, N /7(1 + 2;) ,3 5z - 11 , 2 
j In^ ( 1 - ^) + [-^^^ In' z + In^ z + 



5^-11 



6 

{l + z 
3 

2(1-^: 
3 

(1+2; 



^(2)ln;^ + il±^Li2(l-^) 



12 

(1 + ^) 



■C(2) - 
-C(2) - 



6 
72 



ln(l -z)- 



72 

(1 + ^) 
3 

(1 + ^) 
1152 



43 -5z 
144 



In z 



Si,2(l -z)- ^2^Li3(l - z) 



Wz + 



23-9;2 
1728 



In'^^ 



33 + 5^ h{l + z) 



Li2(l - ^)) In^ z + (ii^Li2(l - z) 



576 



144 



S„(i - .) + 1^ + ?(l±ilc(3) - ^) 1. . + (^5a_l£) 



576 



(i±i)c(2))Li,(l - .) + ii^S,,(l - .) - (l±iiLi3(l - .) 

2(1 + ;2) 



36 

T- /I ^ (1 + ^)q n ^ 5(1 + 2;) 
Li4(l - z) — 82,2(1 - z) — — 81,3(1 - z) 



72 



^(l-) + (i^C(2)-i(i^C(3)|. 



288 



6 



9 



(22) 
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= 32(lnA + ^y-320C(2)(lnA + ^)' + 640C(3)(lnA + ^)' 

+ 240C(4)(^lnA + ^) + 768C(5) - 640C(2)C(3). (23) 

3 Conclusions 



Table 1: Integral I{x) in different approximations. 



X 


Oia) 


0{a') 


0(a3) 


0{a') 


0(a5) 


exponent. 


0.01 


0.99972722 


0.99970278 


0.99970216 


0.99970216 


0.99970216 


0.99970216 


0.1 


0.99713067 


0.99696386 


0.99696237 


0.99696241 


0.99696241 


0.99696241 


0.5 


0.97933805 


0.97871205 


0.97873132 


0.97873149 


0.97873148 


0.97873148 


0.9 


0.91043156 


0.91243286 


0.91253629 


0.91252777 


0.91252803 


0.91252802 


0.99 


0.79019566 


0.80982514 


0.80884689 


0.80886291 


0.80886436 


0.80886422 


0.999 


0.66569598 


0.71915664 


0.71380222 


0.71416672 


0.71415030 


0.71415065 



The expressions obtained for the fourth and fifth order splitting functions were checked to 
satisfy the condition (0). In this way we see the agreement between the iteration procedure 
for Pq^^ (z) and the expansion of the known solution for Our result does also coincide with 
the corresponding expansion of the exponentiated solution from Ref. 0. So, we reproduced 
the known result, but in a different approach; and the higher order splitting functions are 
given explicitly. The exponentiation and order-by-order calculations are complementary to 
each other. As could be seen from Table 1, the numerical difference between the exponentiated 
and non-exponentiated results is negligible, and one may choose safely the approach, which he 
likes. The results of our paper can be used to estimate higher order radiative correction and 
to analyse the numerical difference between exponentiated and order-by-order calculations. 
Here we should note, that in a realistic situation in order to provide a high theoretical precision 
one should take into account also sub-leading radiative corrections, which can be obtained 
only by direct perturbative calculations. 

In Table 1 the values of integral 

1 

I{x) = JdzV!;'{z,Q') (24) 

X 

of the pure photonic part of the non-singlet SF is given for different order approximations for 
Q"^ = 10^ GeV^, L ^ 24.37. For the last column the exponentiated result is obtained by using 
formula (11) from Ref. In Table 2 we present the corresponding values of integrals of the 
splitting functions themselves: 

1 

jW(a;) = J d^p(")(z). (25) 

X 

In the last line of Table 2 we put also the values of the corresponding A-parts. Note, that in 
reality, to simulate the limit A 0, and so to eliminate the dependence on A of numbers for 
integrals I{x) and J^^\x), we used A = 10~^°. 
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Table 2: Integrals of splitting functions J*^")(x) and . 





Oia) 


0{c^) 




0{a') 


0{a^) 


X 




0.01 


-0.0101 


-0.0664 


-0.1854 


-0.1409 


0.5232 


U.i 


-0.1057 


-0.4529 


-0.4487 


1.6579 


2.7010 


n K 
u.o 


-0.7613 


-1.6997 


5.7829 


7.5130 


-99.1465 


0.9 


-3.3002 


5.4338 


31.0395 


-376.7425 


2102.5492 


0.99 


-7.7303 


53.2969 


-293.5939 


708.7696 


11810.0253 


0.999 


-12.3175 


145.1533 


-1606.9698 


16122.5016 


-133803.0981 


A 


p(") 


0.001 


-12.3155 


145.0921 


-1605.5843 


16095.0956 


-133303.9430 



Appendix A 



The Riemann (^-functions are defined as usual: 

= Et^T' C(2) = -, C(3) ^ 1.20205690315959, 

fc=i ^ ^ 

C(4) = ^, C(5) ^ 1.03692775514337 . (A.l) 

Here we define the polilogarithms, which enter into our formulae. We follow the notations 
of Ref . 0, . The general Nielsen's polilogarithm is 



S„,„(z) = , / — ln"-^(x) ln'"(l - xz), (A.2) 

[n — l)!m! J x 







in particular 

Li2(^) = Si,i(^) = - / — In(l-xz), Si,2(^) = ^ / — ln2(l 

J X I J X 



113(2;) = 82,1(2;)= / — ln(x) ln(l - x^;) = / — Li2 

J X J X 


Li4(z) = S3,i(;.) = -i/^ln2(x)ln(l-x^), = J 



1 } dx , r, , „ . . 1 } dx . r..^ . 

/ — In^l-xz) 
b J X 





1 r dx 



S2,2(^) = / — ln(x)ln2(l-xz). (A.3) 
2 J X 







In order to get expansion (|T7D it is convenient to use the following representation for 
F-function: 

1 r^^r.z 



ze 



r(^) „=iV n 
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